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ABSTRACT 


Brouwer and Brouwer-Lyddanes’ use of the Von Zeipel-Delaunay method is employed to 
develop an efficient analytical orbit theory suitable for micro-computers. A succinctly simple 
pseudo-phenomenologically conceptualized algorithm is introduced which accurately and 
economically synthesizes modeling of Drag effects. The method epitomizes and manifests 
effortless efficient computer mechanization. Simulated (Space Telescope) trajectory data is 
employed to illustrate the theory’s ability to accurately accommodate oblateness and Drag 
effects for micro-computer ground based or on-board predicted orbital representation. Real 
(SMM — Solar Maximum Mission) tracking data is used to demonstrate that the theory’s orbit 
determination and orbit prediction capabilities are favorably adaptable to and are comparable 
with results obtained utilizing complex “Definitive Cowell Method” solutions on satellites expe- 
riencing significant Drag effects. 
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AN ECONOMICAL SEMI-ANALYTICAL ORBIT THEORY FOR 
RETARDED SATELLITE MOTION ABOUT AN OBLATE PLANET 


BY: ROBERT A. GORDON 
(NASA/GSFC) 


INTRODUCTION : 

Brouwer* derived a first-order perturbation solution expressing the secular, short 

and long periodic variations in the motion of an artificial satellite about an 

oblate planet. Brouwer obtained separation of all the periodic terms by adapting 
2 

Von Zeipel's technique to modify Delaunay's method for calculating the coeffi- 
cients of the periodic terms through a succession of canonical transformations. 
Delaunay's variables were introduced in order to simplify the canonical expres- 
sions for the equations of motion. Brouwer developed the periodic terms to0GJj} 
and obtained the secular variations to0(^). The resultant formulas are piece- 
wise continuous with singularities existing for certain values of the eccentricity 
and inclination which occur as poles in the algebraic expressions. Thus, the 
equations are valid, except in the regions for which *6 Jet- , and i-£co$*i+ 

sO® l.e., i f/ 2.(,5‘H3 0 , the critical inclination. Lyddane 3 introduced Poincare's 
variables and reformulated Brouwer's expressions as to remove the poles, and thus 
the singularities arising from small eccentricities or inclinations in the Brouwer 
theory. 

This paper is the fruition of an effort to provide an optimal on-board ephemeris 
representation employing an efficient analytical orbit theory suitable for micro- 
computers. Brouwer/Brouwer-Lyddane's method is modified to develop an economical 
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analytical orbit theory for satellite motion about an oblate planet which accommo- 
dates J 3 an£ * parts of the zonal effects with the true argument of latitude 
as the fast variable. The theory is applicable to circular and non-circular satel- 
lites but singular for 2=0 . This is satisfactory for the vast majority of our 
satellite support. The choice of the true argument of latitude as the fast varia- 
ble in difference to Brouwer-Lyddane 1 s choice of the true longitude is a major 
contribution to the economical variation presented here; it simplifies the com- 
putation of the osculating inclination. and portions of the zonal effects 
are considered in the theory in relation to their primary effects on the radial 
and cross-track errors respectfully, and truncated in accordance with economical 

computational consideration. Lyddane remarks that l and 0 must be used for com- 

f t / 4 

and Y in his version; however, as demonstrated by Gordon et al., this 

results in a relative large radial error with respect to Brouwer for moderate 
values of the eccentricity. This can be avoided by evaluating f and X' with £,0* 

h §i 

for moderate values of the eccentricity and with£, 0 for relative low values 

f t 1 

and jC wi th the 

long-period contribution (J 3 ) to the eccentricity. For some orbital parameters, 
this can result in a significant improvement in accounting for intrack error due 
to the oblateness perturbation and compares favorably with respect to the Brouwer 
and Brouwer-Lyddane' s orbit theories for the satellite cases presented in Refer- 
ence 4. The theory presented here is further modified to incorporate a "cheap" 
algorithm which accounts for drag effects semi-analytical ly. A succinctly simple 
pseudo-phenomenologically conceptualized algorithm is introduced which accurately 
and economically synthesizes modeling of drag effects. The method epitomizes and 
manifests effortless efficient computer mechanization. Simulated (Space Tele- 
scope) trajectory data is employed to illustrate the theory's ability to accurately 
accommodate oblateness and drag effects for microcomputer ground-based or on-board 
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predicted orbital representation. Real (SMM-Solar Maximum Mission) tracking data 
is used to demonstrate that the theory's orbit determination and orbit prediction 
capabilities are favorably adaptable to and comparable with results obtained 
utilizing complex "Definitive Cowell Method" solutions on satellites experienc- 
ing significant drag effects. 

FORMULAS FOR COMPUTATION : 

A computational flow diagram of a subroutine with a description of input and out- 
put parameters for that part of the modified Brouwer theory which accounts for 
the oblateness effects is presented in the Appendix. Henceforth, this analytic 
part of the current orbit theory will be represented by the symbolic function 
Bg(t), where t designates the time of theory evaluation. 

The theory is adapted to accommodate retarded motion due to drag by a pseudo- 
physical secular relationship to describe decay in the semi-major axis. This 

representation is inferred phenomenologically from the signature of the semi- 

5 

major axis Locus defined by osculating to mean conversions of state vectors 
of a drag perturbed satellite ephemeris. 

OSCULATING-TO-MEAN CONVERSION : 

5 

Walter's algorithm for osculating to mean conversion is unstable for low e in 
Keplerian space; the apparent instability of the iterative osculating- to-mean 
element conversion is removed by translating the iteration from mean Keplerian 
space to mean Cartesian space. 

Define : 

=• ^ l ) -- Mean Keplerian Elements 

Osculating Keplerian Elements 
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Xs(x',a 

Ys ~ 


Mean Cartesian State Elements 
Osculating Cartesian State Elements 


Given an osculating Cartesian state Y* we determine 


Where j^g represents the Keplerian state two-body functional relationship to 
the Cartesian state. Then employing the iterative algorithm, 

Q**- B 3 (a t*o, Q cl> ) 


or 


LOO 


D ° E. (x'^ 0 ) 

For j * 0, 1, 2, ...» 10 or until the following criterion is satisfied: 

Where £ is some preassigned small positive number. Let this algorithm be 
represented by the symbolic functional relationship. 


Q+-0CY) 


SEMI-MAJOR AXIS DECAY RATE: 


Applying the osculating to mean conversions at one period (f) intervals, we deter- 
mine the semi-major axis decay over M periods, i.e., with 

< «- 0 ,( 3 ) 

Given for i = 1, 2, . . ., M; we compete the mean semi-major axis decay rate 
as 


• « 




•f n o 

Pit *" 

n 
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ORBIT PROPAGATION WITH "Bq" : 

To update the orbital elements to time (at with the Bg theory, we assume 
the orbital elements remain constant over one period and rectify the theory's con- 
stants at one period intervals (with the a decay rate) up to the Nth period where 

N = id • 

Employing the above iterative method we have 

i = o-i-. + «-P H 
(f°)i ~ 

Evaluating the secular part of Bg we obtain 

j m N 

Then at timedithe osculating elements are given by evaluating the full Bg 
theo»7 with 

8, (pa), At -T^) 

With T^ = N x P. Let us represent the semi-analytic theory with the recti fi- 
cation algorithm for retarded motion symbolically by "Bg." 


TJ = 



U. 

of* 


Where 


TRAJECTORY DATA : 

Trajectory data, i.e., osculating state vectors are used in a simulation to demon- 
strate the Bg orbit theory capability in representing retarded satellite motion 
about an oblate planet. It has been proposed that a secular analytic orbit theory 
be employed for the on-board ephemeris representation of the Space Telescope. The 
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Space Telescope can experience significant drag effects over a three-day span. 
State vectors generated from a sample set of Space Telescope elements demon- 
strates the Bg theory's superior ability to represent the Space Telescope 
ephemeris. This will be demonstrated by three steps in the simulation. 


Step No. 1 : 


A comparison is made between two Cowell ephemeris generations at two-hour inter 
vals, the "Truth Ephemeris" with drag included with a drag model constant C n = 


2.0 versus the Cowell ephemeris with C| 
a three-day span. The Space Telescope 

CL = 6778.140 km 

e = o.ooi 

1 = 28.2 degrees 
tl = 19.78 degrees 
3=0 degree 
' 1=0 degree 


( = 0, i.e., no drag consideration over 
epoch elements is defined as: 

A = 117.6 m 2 
to = 10134 kg 

= 393.222 km 
>| a = 406.778 km 


Note: Table No. 1 — the maximum total error growth realized was (262.88 km). 


Step No. 2 : 

An analytical method without drag model effects is fitted to the "Truth Ephemeris" 
over a three-day span with the "Truth Ephemeris" state vectors as observation 
data at two-hour intervals for a Differential Correction of the epoch mean ele- 
ments of the analytical orbit theory. The analytical theory used is the Brouwer- 
Lyddane which includes periodic terms. 

Note: Table No. 2 — the post trajectory data DC compare at the two-hour fre- 
quency yields a maximum total error of (43.71 km). 
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Step No. 3 : 


Step No. 2 was repeated employing the Bg orbit theory. 

Note: Table No. 3 — the post trajectory data DC compare at the two-hour fre- 
quency yields a maximum total error of (1.30 km). 

TRACKING DATA : 

Real tracking data demonstrates the Bg orbit theory's favorable orbit determina- 
tion and prediction capabilities. An orbit determination for a number of differ- 
ent epochs employs real (SMM-Solar Maximum Mission) tracking data over a two-day 
span to differentially correct the epoch state and drag model parameter for a 
"Definitive Cowell Method" and the epoch mean elements for the Bg orbit theory. 
The predicted ephemeris of the "Definitive Cowell Method" and the Bg orbit theory 
is then compared with a series of state solutions determined over a two-day DC 

arc at two-day intervals. A table of the comparable response of the Cowell and 

• 

Bg method is presented in Table 4. 

CONCLUSIONS : 

1. The osculating to mean algorithm described herein provides an accurate 
first-order estimate to the semi -major axis decay rate. 

• 

2. Tables 1 and 2 graphically demonstrates that the Bg orbit theory can 
accurately accommodate significant drag effects on orbital motion and 
that the a parameter can absorb virtually all of the significant drag 
effects. 

3. Table 4 implies that the Bg theory is competitive with a "Definitive 
Cowell Method" in a least squares batch filter orbit determination for 
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satellites experiencing significant drag effects. Thus this can expand 
the class of satellites which can be operationally supported with a semi- 
analytic orbit theory. 

• 

4. The Bg theory is suited for a ground-based or on-board microcomputer appl 
cations, providing an orbital ephemeris generation which does not require 
a density table or analytic density model. 

RECOMMENDATIONS : 

1. Some adaptation of the Bg theory should be employed for the on-board 
ephemeris representation of the Space Telescope. 

2. Develop the state transition matrix for a truncated secular version of 
Bg for Karman Filter state estimation applications on microcomputers. 

3. Bg be adapted by those various sites who require orbit ephemeris genera- 
tion but does not have an orbit determination capacity. The mean orbital 
constants and the a parameter can be determined and distributed by 
Goddard for satellites of interest as are now the Brouwer Mean Orbital 
elements. This would lead to a uniform method at each of the various 
sites who require such elements with a significant improvement in ephem- 
eris representation for satellites experiencing significant drag effects. 
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POSITION DIFFERENCES 
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POSITION DIFFERENCES 

POSITION DIFFERENCES 

OF 


(METERS) 



DEGREES 

SEC.S 

(METERS) 


(KM. ) 

(METERS) 

(OEG. ) 

DATA 

• * * “ 

* * * * 

PPM 

( 10** 

3) 

‘ - — - - 

• • * * 

“***“■ 

TOTAL 

• • * * 

IN 

ALONG 

CROSS 

TRUE 

YYMMDD 

HHMMSS 

DA 

OE 

DI 

DELH 

DELG 

DELM 

DELTIME 

ERROR 

RADIAL 

TRACK 

TRACK 

TRACK 

ANOMALLY 

800101 

0 . 

0.00 

0.0 

0.00 

-0.00 

-0.00 

0.00 

-0.0000 

0.00 

0.00 

-0.00 

-0.00 

0.00 

360.00 

800101 

20000. 

83.59 

24.41 

-1.53 

-2.27 

-0.05 

0.06 

0.0102 

240.24 

-80.42 

0.08 

78. 15 

-212.47 

350.43 

800101 

40000. 

-127.00 

-8.17 

-0.96 

-3.62 

-0.94 

0.94 

0.0927 

721.76 

-69.67 

0.71 

711.72 

-97.69 

1.24 

800101 

60000. 

17.94 

-10. 19 

-0. 11 
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0.57 

-0.56 

0. 1455 

1 132.85 
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1 . 12 
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800101 

80000. 

-178. 12 

33.25 
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-1.98 
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2410. 19 
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4.99 

4991 . 17 
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800101 

120000. 

-143.57 

-62.50 

-0.59 

0.54 

-0.57 

0.62 

0.7995 

6143.41 

276.48 

6. 14 

6136.94 

54.85 

19. 13 

800101 

140000. 

-267.31 

30.37 

-1.96 

0.80 

3.47 

-3.40 

1 .01 15 

7776.03 

-375.61 

7.76 

7764.87 

-179.85 

344.47 

800101 

160000. 

-208.34 

54.51 

-1.34 

-3.09 

-1.15 

1.25 

1.5163 

11659.37 

-565.02 

11.64 

1 1643.32 

-234.01 

1.41 

800101 

180000. 

-300. 10 

-18.26 

0.49 

-3.84 

-2.48 

2.60 

1 . 9604 

15050.25 

-134.43 

15.05 

15049.23 

-112.15 

7.71 

800101 

200000. 

-328.03 

-16.43 

1.80 

-2.37 

1 . 17 

358.99 

2.4166 

18554.05 

-190.65 

18.55 

18552.48 

147.85 

359.38 

800101 

220000. 

-463.20 

10.98 

1.76 

-0.07 

1.43 

-1.23 

3.0486 

23409.93 

-489.38 

23.40 

23404.04 

191.29 

355.61 

800102 

0 . 

-346.01 

9. 18 

-0.61 

-0.56 

-1.21 

1.45 

3.7498 

28793.71 

-355.66 

28.79 

28791.50 

22.34 

357.94 

800102 

20000. 

-387.82 

-30. 1 1 

-0.96 

-5.09 

0.87 

-0.58 

4.3352 

33282.49 

-116.42 

33.28 

33280.88 

306.77 

9.53 

800102 

40000. 

-535.87 

23.41 

-1.00 

-6.99 

2.68 

-2.35 

4.9812 

38241.07 

-543.28 

38.24 

38237.20 

-24.84 

346.46 

800102 

60000. 

-450.74 

25.05 

-0.78 

-4.87 

-1.39 

1.77 

5.8524 

44943.41 

-477. 19 

44.94 

44940.08 

-267.50 

359.53 

800102 

80000. 

-558. 16 

-59.85 

-0. 11 

-2.89 

-0.35 

0.78 

6.5768 

50486.37 

41.35 

50.49 

50486.29 

78.00 

12.65 

800102 

100000. 

-505.33 

55.42 

-1.28 

-1.15 

4.61 

-4.12 

7.3302 

56277.95 

-564.37 

56.28 

56275. 10 

-52.92 

348.25 

800102 

120000. 

-582. 19 

70.86 

-1.36 

-1.87 

-2. 18 

2.72 

8.4557 

64931.37 

-728.67 

64.93 

64927.05 

-171.64 

4.62 

800102 

140000. 

-621.43 

-64.72 

-1.83 

-3.82 

356.46 

-355.86 

9.3871 

72060.31 

207.04 

72.06 

72059.81 

169.50 

3.15 

800102 

160000. 

-672.69 

-34.59 

-1.31 

-3.82 

4.05 

356.62 

10. 1744 

78116. 10 

22.44 

78. 12 

781 15.85 

196.97 

357. 14 

800102 

180000. 

-646.32 

61.29 

1.11 

-5.38 

1.95 

-1.21 

11.3165 

86879.22 

-509. 12 

86.88 

86877.73 

-12.47 

2.72 

800102 

200000. 

-739.65 

13.62 

1.48 

-3.37 

-1.51 

2.33 

12.6381 

97034.61 

- 166.98 

97.03 

97034. 12 

-259.67 

354.53 

800102 

220000. 

-692.36 

-25.74 

1 . 37 

-2.44 

0.91 

-0.01 

13.8261 

106155.76 

313.94 

106. 16 

106155.28 

46.86 

6.65 

800103 

0 . 

-734.94 

31.32 

-1.29 

-5.37 

3.78 

-2.79 

15.0265 

1 15351 .96 

48.73 

115.35 

1 15351 .91 

-95.41 

355.39 

800103 

20000. 

-935.52 

19.92 

-1.56 

-5.78 

-0.64 

1.71 

16.4412 

126243.38 

95.53 

126.24 

126242.99 

-296.90 

357.94 

800103 

40000. 

-715.54 

-34.64 

0.02 

-8.97 

-1.16 

2.32 

17.7319 

136131.57 

914.00 

136.13 

136128.31 

230.72 

8.14 

800103 

60000. 

-1029.21 

-3.00 

-0.89 

-7.32 

5.57 

-4.32 

18.9784 

145708.60 

606. 19 

145.71 

145707. 18 

216.61 

352.29 

800103 

80000. 

-868.53 

73.62 

0.40 

-5.34 

-0.52 

1 .86 

20.6096 

158243.87 

505.38 

158.24 

158242.86 

-252. 15 

5.82 

800103 

lOOOOO. 

-945.72 

-78.58 

-0.56 

-6.83 

-3.02 

4.46 

22. 1222 

169827.86 

1682.91 

169.82 

169819.52 

12.59 

357.00 

800103 

120000. 

-1120.45 

-48.86 

-1.89 

-1.87 

6.08 

355.45 

23.4078 

179744.87 

1619. 19 

179.74 

179737.56 

63.59 

355.90 

800103 

140000. 

-870.93 

108.63 

-1.23 • 

-8.39 

4.08 

-2.44 

25. 1243 

192892.69 

1114.23 

192.89 

192889. 14 

-357.96 

8.79 

800103 

160000. 

-1203.34 

-29.79 

-1.36 

-9.34 

-3.35 

5. 10 

26.9431 

206860.50 

2062.75 

206.85 

206850.07 

-242.30 

352.43 

800103 

180000. 

-lOOO. 10 

-55.88 

2.04 

-7.64 

2.47 

-0.61 

28.5222 

219028.26 

2918. 1 1 

219.01 

219008.62 

296.97 

2.94 

800103 

200000. 

-1209.92 

11.27 

1.93 

-9.26 

7.35 

-5.36 

30.2842 

232502.55 

2673.50 

232.49 

232487.01 

281.67 

2. 13 

800103 

220000. 

-1187.26 

48.39 

-0.02 

-3.44 

0. 19 

1 .92 

32.3532 

248429.29 

3026.44 

248.41 

248410.78 

-190.99 

356.58 

800104 

0 . 

-1008.66 

-31.85 

0.08 

-11.63 

-0. 1 1 

2.34 

34.2400 

262918.02 

4321.32 

262.88 

262882.30 

329.56 

2.74 


Table No. 1 -- Cowell Drag Versus Cowell No Drag 
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-0.52 

-3.00 

-2.03 

1.86 

-2.5524 

19592.75 

-19.64 

-19.59 
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327.99 

-21.08 

-21081.00 

-105.67 

359.53 

800102 

80000. 

-1.57 

-45.45 
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Table No. 2 -- Cowell Drag Versus Brouwer-Lyddane 
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